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Abstract . 

"^Transient  acoustic  wave  propagation  is  analyzed  for  the  case  of 
an  unlimited  plane-stratified  fluid  having  constant  density  and  sound 
speed  c(y)  at  depth  y  given  by  the  Epstein  profile 


£ 


c~^(y);  *  K  sech2  (y/H)  +  L  tanh  (y/H)  +  M 


.  Cj 


The  acoustic  potential  is  a  solution  of  the  wave  equation 

Qt  *  it*  4  •  *4  tiff'd  ~.v  -  •  .  r  ■> 

Diju  ~  c2  (y)  (of)u  +  D|u  +  D*u)  -  f (t,x,y) 


f  J  f  ^  *  -  .  .  f 


where  x  »  (xj,x2)  are  horizontal  coordinates  and  f(t,x,y)  characterizes 
the  wave  sources.  The  principal  results  of  the  analysis  show  that  u  is 
the  sum  of  a  free  component,  which  behaves  like  a  diverging  spherical 
wave  for  large  t,  and  a  guided  component  which  is  approximately  localized 
in  a  region  ] y  —  yQ |  <  h  and  propagates  outward  in  horizontal  planes 
like  a  diverging  cylindrical  wave. 


C  t 


§1.  Introduction. 


This  paper  presents  an  analysis  of  transient  acoustic  wave 
propagation  in  a  stationary  unlimited  plane  stratified  fluid  with 
constant  density  and  sound  speed  c(y)  at  depth  y  defined  by  the  Epstein 
profile 

c(y)  ■  [K  sech2  (y/H)  +  L  tanh  (y/H)  +  M]-1^2  (1.1) 

where  H,  K,  L  and  M  are  constants  such  that  c(y)  >  0.  The  acoustic 
field  is  characterized  by  a  real-valued  potential  u(t,x,y)  that  satis¬ 
fies  the  wave  equation  [8] 

D2u  -  c2 (y) (D2u  +  D2u  +  D2u)  -  f(t,x,y)  (1.2) 

where  t  is  a  time  coordinate,  x  »  (Xj,x2)  are  Cartesian  coordinates  in  a 
horizontal  plane,  f(t,x,y)  is  a  function  that  characterizes  the  wave 
sources  and  »  3/3xj,  Dt  ■  3/3t,  -  3/3y. 

The  sound  speed  profiles  (1.1)  were  introduced  by  P.  S.  Epstein 
[2]  who  discovered  that  (1.2)  with  this  choice  of  c(y)  can  be  integrated 
by  means  of  hypergeometric  functions.  This  fact  is  of  interest  in 
theoretical  acoustics  because  the  Epstein  profile  provides  an  example  of 
the  physically  interesting  phenomenon  of  an  acoustic  duct  when  c(y)  has 
a  minimum.  This  case  is  characterized  by  the  parameter  values  [3,4] 

H  >  0,  M  >  L  >  0,  K  >  L/2  (1.3) 

—  1/2 

For  these  values  the  limits  c(±°°)  -  (M  ±  L)  are  finite  and  c(y)  has 
a  unique  minimum  c(y0)  <  c(°°)  <  c(-°°)  at  y0  »  H  tanh-1 (L/2K) . 
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The  integration  of  equation  (1.2)  below  is  based  on  the  spectral 
theory  of  the  Epstein  operator  A  in  the  Hilbert  space  X  ■  L2 (R3 ,c-2 (y)dxdy) 
as  developed  in  [3,4]  and  the  notation  and  results  of  that  paper  are  used. 
Equation  (1.2)  is  interpreted  as  the  equation 

D|u  +  Au  -  f  (t, •) ,  t  £  R,  (1.4) 

for  a  function  t  -*■  u(t,*)  £  X.  The  wave  sources  are  assumed  to  act 
during  the  time  interval  [0,T],  so  that  supp  f  c  [0,T] .  The  correspond¬ 
ing  acoustic  wave  is  the  solution  of  (1.4)  that  satisfies  the  initial 
condition 


u(t,*)  -  0  for  all  t  <  0  (1.5) 

The  solution  is  given  by  Duhamel's  integral 


u(t,*) 


{A-1/2  sin  (t  -  x)  A 1/2  }  f  (t ,  •  )dx ,  t  >  0 
o 


(1.6) 


Indeed,  if  f  €  C([0,T],-K)  then  (1.6)  is  the  unique  "solution  with  finite 
energy"  of  [9],  while  if  f  e  C([0,T],  D(Al/2  ))  then  (1.6)  is  the  "strict 
solution  with  finite  energy."  In  addition,  if  f  £  C([0,T],  D(A-1^2  )) 
then 


u(t,x,y)  -  Re  (v(t,x,y)}  (1.7) 

where  v(t,*)  is  the  complex-valued  potential  defined  by 


t 

v(t,*)  -  i  exp  {-it  A1/2  }  A~1/2  exp  {it  Al/2  }  f(x,*)dx 


In  particular. 


3 


v(t,*)  a  exp  {-it  A1^2  }h,  t  >  T 


(1.8) 


where 


h  -  i  A'1/2 


fT 

exp  {it  A1^2  }  f(x,*)dT 
o 


(1.9) 


The  initial-value  problem 


D2u  +  Au  *  0  for  all  t  >  0  (1.10) 

u(0)  -  f,  Dtu(0)  -  g  (1.11) 

can  be  treated  by  the  same  formalism.  Indeed,  if  f  6  D(A*/2  )  and 
g  e  D(A-1/2  )  then  the  solution  of  (1.10),  (1.11)  is  given  by  (1.7), 

(1.8)  where  h  -  f  +  i  A*l/2  g  €  D(A1/2  )  (cf.  [11,  Ch.3]). 

The  integral 

E(u,K,t)  -  f  { (DjU) 2  +  (D2u)2  +  (Du)2  +  c-2 (y) (D  u)2}dxdy  (1.12) 
jk  y 

may  be  interpreted  as  the  energy  of  the  acoustic  field  u  in  the  set 
K  c  R3  at  time  t.  Moreover,  A  is  the  self adjoint  operator  associated 
with  the  Dirichlet  integral 


A(u,v) 


2 

I 

j-1 


(R1) 


®y-y)L, 


(R’> 


in  the  sense  of  T.  Kato's  theory  of  sesquilinear  forms  [5].  Indeed,  if 
the  domain  of  A  is  D(A)  *  L2(R3),  the  first  Sobolev  space,  then  A  is 
precisely  the  associated  operator  of  Kato's  theory.  It  follows  from 
Kato's  second  representation  theorem  [5,  p.  331]  that 


D(A1/2  )  -  l4(R3) 


(1.13) 


4 


and 


|AI/2  u|£.  -  A(u,u) 


2 

I 

j-1 


'V 


+ 


■v 


l2(r3) 


(1.14) 


Hence  the  total  energy  satisfies 

E(u,R3 , t)  -  SA1/2  ul2.  +  JD^I2.  (1.15) 

If  h  €  D(A3^2  )  and  u  is  defined  by  (1.7),  (1.8)  then  a  simple  calculation 
shows  that 

E(u,R3  ,t)  =  llAl/2  h!2.  for  all  t  >  T  (1.16) 

The  analysis  of  the  structure  of  the  acoustic  potential  (1.7), 
(1.8)  presented  below  is  based  on  the  eigenfunction  expansion  of  [4].  The 
orthogonal  projections  in  3C 


Pf  »  P+  +  P_  +  P0  (1.17) 

and 

00 

p  «  z  Pk  (1.18) 

8  k-1 

defined  by  the  eigenfunction  expansion  provide  a  decomposition 

vi(t,*)  -  uf(t,*)  +  ug(t,*) 

into  orthogonal  partial  waves 

uf(t,*)  -  Pf  u(t,*),  ug(t,‘)  *  Pg  u(c»‘) 

The  first,  called  the  free  component,  will  be  shown  to  behave  for  large 
times  like  a  diverging  spherical  wave  in  a  homogeneous  fluid.  The 
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second,  called  the  guided  component,  will  be  shown  to  be  approximately 
localized  near  the  plane  y  ■  y0  and  to  propagate  outward  in  horizontal 
planes  like  a  diverging  cylindrical  wave.  This  second  component  shows 
the  profound  effect  of  an  acoustic  duct  on  transient  acoustic  waves.  It 
is  absent  when  c(y)  has  no  minimum. 

Transient  wave  propagation  in  the  analogous  but  simpler  case  of 
the  Pekeris  duct  (c(y)  piece-wise  constant)  was  analyzed  by  Wilcox  in 
[12]  and  [13].  This  paper  is  a  sequel  to  the  report  [10]  of  Wilcox  and 
the  article  [4]  of  Guillot  and  Wilcox.  In  [10]  the  special  case  of  a 
symmetric  Epstein  profile  (L  =  0)  was  analyzed  without  detailed  proofs. 
[4]  presented  a  complete  spectral  analysis  of  the  Epstein  operator.  Here 
the  results  of  [4]  are  used  to  treat  the  general  Epstein  duct  and  to 
supply  the  proofs  that  were  twitted  in  [10].  Some  of  the  results  of 
this  paper  were  announced  in  [1]. 

The  remainder  of  this  paper  is  organized  as  follows.  In  §2  the 
eigenfunction  expansion  of  [4]  is  reformulated  to  provide  a  convenient 
starting  point  for  the  analysis  of  the  free  component  uf.  The  behavior 
of  u^(t,*)  for  t  *►  “  is  calculated  in  §3.  The  justification  of  these 
calculations  is  technically  more  difficult  than  for  the  Pekeris  profile 
treated  in  [13]  and  proofs  of  the  results  are  presented  here  for  the 
first  time.  §4  presents  a  calculation  of  the  asymptotic  behavior  for 
t  -*■  00  of  the  guided  component  Ug(t,*).  55  presents  applications  of  the 

results  of  §53  and  4  to  calculating  asymptotic  distributions  of  energy  for 
large  times.  The  proofs  of  the  results  in  §54  and  5  are  essentially  the 
same  as  those  for  the  Pekeris  profile,  given  in  [13],  and  are  therefore 
omitted.  A  version  of  the  method  of  stationary  phase  for  oscillatory 
integrals  containing  parameters  in  given  in  an  Appendix. 


The  generalized  eigenfunctions  ^+(x,y,p,X)  and  tj>0  (x,y,p,A)  were 
defined  in  [4]  as  multiples  of  certain  hypergeometric  functions.  The 
normalizing  constants  a+(p,A)  and  a0(p,A)  are  determined  only  up  to 
factors  of  modulus  1.  In  [4]  these  factors  were  chosen  to  make  a+(p,A) 
and  a0(p,A)  real  and  positive  [4,  pp.  92-93],  Here  it  will  be  convenient 
to  renormalize  to  make  the  constants  c+(p,A)  and  c0(p,A)  of  the 
asymptotic  forms  [4,  (1.7),  (1.8),  (1.11)]  real  and  positive.  Calcula¬ 
tion  of  c+(p,A)  and  c0(p,A)  when  a+(p,A)  and  a0(p,A)  are  defined  as  in 
[4]  shows  that  the  quantities  (4TTq+)1^2  c+(p,A)  *  exp  {i  ¥+(p,A)}  and 
(4irq+)1^2  c„(p,A)  =  exp  {i  4'„(p,X)}  have  modulus  1.  Hence,  if  and 
are  renormalized  by  replacing  a+(p,A)  and  aQ(p,X)  by 
a±(p,A)  exp  {-i  'F±  (p,X) }  and  a0(p,A)  exp  {-i  'P<)  (p »X) }  then 

(4uq±)l/2  c±(p,X)  -  1,  (4irq+) c0 (p,X)  ■  1  (2.1) 

Explicit  expressions  for  exp  {i  V+(.p,\)}  and  exp  {i  4'0(p,X)}  as 
quotients  of  products  of  T-functions  are  easily  obtained  but  will  not  be 

needed  here.  The  normalization  (2.1)  is  employed  in  the  remainder  of 

this  paper. 

The  generalized  eigenfunctions  \p+  and  ip0,  renormalized  as  above, 

satisfy 

^(y.P.A)  -  (4TTq+)-1/z  T±(y,p,X)  exp  { +i  q_j-}  (2.2) 

^0(y»P,X)  =*  (4TTq+)_l/z  T0  (y,p,X)  exp  {0  y/H}  (2.3) 

where 

i-HLi- JLui Ci  bn^'JL  K«“i*QT  i  1 
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T+(y,p,A)  -  (47^) 1/2  a+(p,A)  exp  {-if^p.A)  +  iq_y)  4>3(y,p,A) 

T_(y,p,A)  »  (4nq_)1/2  a_(p,A)  exp  {-i'F^p.A)  -  iq+y}  $4(y,P,A) 

T0(y.P»A)  =  (4irq+)1/2  a0(p,A)  exp  {-i'4'0(p,A)  -  8y/H}  4>x  (y ,p ,A) 

and 

11m  T+(y,p, A)  =>  T±(p,A)  (2.4) 

y-+jpo 

lim  T0(y,p,A)  =  T0(p,A)  (2.5) 

oo 

Similarly,  well-known  identities  for  the  hypergeometric  functions 
imply  [3] 

>My>P.A)  =  (4irq+)-1/z  {I  (y,p,A)  exp  (+iq+y)  +  R+(y,p,A)  exp  (±iq+y)} 

(2.6) 

=  (47rq+)_1/2  {l0  (y,P,A)  exp  (-iq,y)  +  RQ(y,p,A)exp  (iq+y)} 

(2.7) 

where 

I+(y>P,*)  *  (2/H) 1/7  exp  (iq+y)  <t>2(y»P»A) 

I_(y»P»^)  ”  (2/H) 1/2  exp  (-iq_y)  4>1(y,p,\) 
l0(.y,pA)  -  (2/H)l/2  exp  (iq+y)  4>2  (y ,  p ,  A) 

R+(y.P»^)  ■  (2/H)1/2exp  (-iq+y)  R+(p,A)  <{>4(y,p,A) 

R_(y»P.^)  =*  (2/H) 1,2  exp  (iq_y)  R_(p,A)  <j>3(y,p,A) 

R0(y.P^)  =*  (2/H)1/zexp  (-iq+y)  R0(p,A)  <f>4(y,p,A) 
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and 


lim 

i+(y.p,A) 

=  1 

(2.8) 

y-^+OO 

lim 

i0(y.p.A) 

=  1 

(2.9) 

y-*+°° 

lim 

R+(y>p»A) 

=  R±(p,A) 

(2.10) 

y-^+OO 

lim 

R0(y>p.A) 

*  R0(P»A) 

(2.11) 

y-*+°° 

calculated  expressions 

for  R±(p,A) , 

R0(p,A),  T+(p,X)  and 

T„(p,A)  as  quotients  of  products  of  T-functions  were  given  in  [3]  and 
will  not  be  needed  here.  They  imply  that 

q+  \K\2  +  q.  |t+|2  *  q+ »  |R„I  *  1  (2.12) 

-  -  +  -  - 

The  construction  of  u^(t,*)  =  Pj  u(t,’)  will  be  simplified  by 
representing  *  P+  +  P_  +  P0  by  means  of  a  single  family  of  generalized 
eigenfunctions,  rather  than  the  three  families  i//+,  ip_  and  i|/0  .  This  may 
be  motivated  by  noting  that  \\t+,  ip_  and  ip0,  collectively,  represent  the 
response  of  the  Epstein  fluid  to  the  incident  plane  waves 
exp  {i(p  •  x  -  q  y) } ,  (p,q)  S  R3.  To  see  this  consider  the  mappings 


(p,q)  -  X+(p,A)  -  (p,q+(p,X)), 

(p,X)  e  a 

(p.q)  3  x0(p»^  *  (p,q+(p.X)), 

(p,A)  g  Qq 

(p , q)  *  X_(p,X)  -  (p , -q_ (p ,X) ) , 

(p,X)  €  Q 

X+  is  an  analytic  transformation  of  SI  onto  the  cone 


10 


c+  =  { (p,q) :  q  >  a|p  I } 

where 

a  =  ((c(-»)/c(«))2  -  1)1/2  >  0 

Similarly,  X0  is  an  analytic  transformation  of  ft0  onto  the  cone 
c0  =  ( (p ,q) :  0  <  q  <  a  I p I } 

and  X_  is  an  analytic  transformation  of  ft  onto  the  cone 

c_  =  { (p,q) :  q  <  0} 

Thus,  the  asymptotic  forms  of  \f>+  and  ip0  for  y  -*■  ±»  [4,  (1.7),  (1.8), 
(1.11)]  show  that  i|)+(x,y,p, A)  with  (p,A)  6  ft  is  the  response  of  the 
medium  to  a  plane  wave  exp  {i(p*  x  -  q  y)}  with  (p,q)  £  C+,  \p0  (x,y,p,\) 
is  the  response  to  a  plane  wave  with  (p,q)  £  C0  and  i]j_(x,y,p,A)  is  the 
response  to  a  plane  wave  with  (p,q)  e  C  .  Note  that 


R3  =  C+  U  C„  U  C_  U  N 
where  N  is  a  Lebesgue  null  set. 

The  interpretation  of  4^,  ip_  and  ip0  given  above  suggests  defining 
the  composite  eigenfunction 


<f>+(x,y,p,q) 


exp  (ip  •  x)  4>+(y>P>q)  ,  (p , q)  e  c+  u  c0  U  C_  (2.13) 


by 


4>+(y,p,q) 


(2.14) 

(2q)  1/2  c(”)  i|J+(y,p,A),  ( p , A )  =  X”1  (p , q)  ,  (p,q)  €  C+ 

=  <  (2q)1/z  c(«)  (y,  p ,  A)  ,  (p,A)  =  X~L(p,q),  (p,q)  €  C0 

(2|  q|  )  1/2  c(-«)  i|/_(y,p,A),  (p,A)  “  X~l  (p,q)  ,  (p,q)6C_ 


The  terms  (2q)  ^2  c(°°)  and  (2  Jq  j )  c(-«°)  are  normalizing  factors.  Note 

that  2q  c2(°°)  is  the  Jacobian  of  X+*  and  Xq1  and  2q  c2(-«)  is  the 
Jacobian  of  X_*. 

The  eigenfunction  expansion  of  [4]  will  be  reformulated  in  terms 
of  <j>+.  To  this  end  let  h  S  3C  and  note  that  h^  ■  P^h  *  P+h  +  P_h  +  P„h 
has  the  representation 

hf (x,y  )  =  (x,y,p,X)  h  (p,X)  dpdX 

r  + 


+  <|j0  (x,y,p,X)  h0(p,X)  dpdX 


+  ip  (x,y,p,X)  h  (p,X)  dpdX 


where  the  integrals  converge  in  JC.  The  JC-lim  notation  will  be  suppressed 
for  brevity.  Changing  the  variables  in  the  three  integrals  by  means  of 
X+,  X0  and  X_,  respectively,  gives 


hf(x,y)  =  <J>+(x,y,p,q)  h+(p,X)  c(“)  (2q)1/2dpdq 

+  <P+(x,y.p,q)  h5(p,X)  c(°°)  (2q)1/zdpdq 

r 


+  4>.  (x.y.p.q)  h  (p,X)  c(-»)  (2 1  q|  )  1/2  dpdq 

JC 


<J>+(x,y,p  ,q)  h  (p,q)  dpdq 
R3 


where 
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X 


*(p,q) 


c2(°°)(|p|2  +  q2), 
c2 (-“) ( | P | 2  +  q2), 


q  >  0 

> 


q  <  0  J 


(2.15) 


and 


h+(p,q) 


< 


(2q)1/2c(°°)  h+(p,X(p,q)) , 
(2q)1/2c(“)  ho  (P»X(p,q)>  , 

<2  |q  |  )  1/2  c(-»)  h_(p,X(p,q)), 


(P,q)  e  c+ 
(p»q)  e  c0  > 
(p,q)  e  c_ 


(2.16) 


It  is  easy  to  verify  by  considering  the  three  cones  C+,  C0  and  C 
separately  that 


h.  (p,q) 


4’4.(x.y.P>q)  h(x,y) 


c  2  (y) 


dxdy 


where  the  integral  converges  in  L2(R3).  Moreover,  it  can  be  shown  by 
direct  calculation,  using  the  Parseval  formula  of  [4],  that 

11  hf  II  ^  -  llh+»L^ 

These  considerations  and  the  eigenfunction  expansion  theory  of  [4,  §5] 
imply 

Theorem  2.1.  For  all  h  S  JC  the  strong  limit 


h,(p,q)  “  L2(R3)-lim  $,(x,y,p,q)*  h(x,y)  c-2 (y)  dxdy  (2.17) 

+  j  x  1 2+y2<M2  + 

3  A 

exists.  Moreover,  the  mapping  3C  L2  (R  )  defined  by  f)+h  ■  h+  is  a 

partial  isometry  such  that 

tt.Q*  -  1  and  n*&,  -  P, 

+  +  +  +  r 


> 


» 


» 


» 


» 


# 


I 


t 
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*  '”N 

The  adjoint  mapping  h+  is  given  by 

b^ (x,y)  =  3C-lim  <f>.  (x,y,p,q)  h  (p,q)  dpdq  (2.18) 

f  m~>  J  |p|2+q25M2  +  + 

Finally,  is  a  spectral  mapping  for  A  in  the  sense  that  for  all  h  e  D(A) 
one  has 


(fi+Ah)  (p , q)  -  X(p,q)  ft+h(p,q) 


(2.19) 


where  X(p,q)  is  defined  by  (2.15). 

Note  that  these  results  are  simply  a  reformation  of  the  results 
of  [4,  §5]  and  not  a  new  theorem. 

It  is  important  for  the  calculations  of  §3  to  recognize  that 
another  family  of  generalized  eigenfunctions  of  A  is  defined  by 


4>_(x,y,p,q) 


<f>+(x,y,-p,q) 


* 


(p » q)  e  r3  -  n 


It  is  clear  that  A<f>_  =  X  (p , q) <{)_  and 

<t>_(x,y,p,q)  =  (2ir)-1  exp  (ip  •  x)  4>_(y,P,q) 


Moreover, 

^.(y.p.q)  “  0+(y»p.q)*  (2. 20) 

because  <j>+(y,p,X),  $0(y,p,X)  and  therefore  <f+(y,p,q),  depend  on  p  through 
|p|  alone.  The  asymptotic  behavior  of  <f>+  and  <f>_  for  y  00  may  be  seen 
from  (2.2),  (2.6),  (2.7),  (2.14)  and  (2.20).  It  is  given  by 


♦ 
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c (“) { exp  (-iqy)  +  R+(p,X)  exp  (iqy)},  (p,q)  £  C+ 

4>+(y.P.q)  ~  (2ir)“1/2  <  c(«){exp  (-iqy)  +  R„(p,X)exp  (iqy)},  (p,q)  £  C0 

c(_oo)  T_(p,X)  exp  (iq+y),  (p,q)  £  C_ 

and 

c (°°) { exp  (iqy)  +  R+(p,X)*exp  (-iqy)},  (p,q)  e  C+ 

<t>_(y.P>q)  -  (2tt)"1/2  <  c(°°){exp  (iqy)  +  R„(p,X)*exp  (-iqy)},  (p,q)  6  C0 

c(-°°)  T_(p,X)*exp  (-iq^y) ,  (p,q)  e  C_ 

These  relations  clearly  imply  that  (y , p , X)  is  not  simply  a  multiple  of 
4>+(y,p,X).  By  contrast,  the  guided  mode  eigenfunctions  have  the  symmetry 
property 

\(y,p)  “  I}»k(y,-p)*»  k  =  1>2,  •  •  •  (2.21) 

because  they  are  real-valued  and  depend  on  p  only  through  | p | . 

The  family  <j)_(x,y,p,q)  ,  (p,q)  £  R3  -  N,  is  a  second  family  of 
generalized  eigenfunctions  for  A  that  spans  the  reducing  subspace 

=  Pj3C.  In  fact,  the  following  exact  counterpart  of  Theorem  2.1  holds. 
Theorem  2.2.  For  all  h  e  JC  the  strong  limit 

h  (p, q)  *  L2(R3)-lim  0  (x,y,p,q)*  h(x,y)  c_2(y)  dxdy 

M-«°  J|x|2+y2<M2 

exists.  Moreover,  the  mapping  X  ■*  L2(R3)  defined  by  ft_h  -  h_  is  a 

partial  isometry  such  that 

•k  k 

£2_&_  ■  1  and 

*A 

The  adjoint  mapping  is  given  by 


IftH  fT"  ' 
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hf(x,y)  =  3f-lim  <f>_(x,y,p,q)  h_(p,q)  dpdq 

M-*»  J  |  p|  2+q2<M2 

Finally,  Q_  is  a  spectral  mapping  for  A  in  the  sense  that  for  all 
h  €  D(A)  one  has 


(fi_Ah)  (p,q)  =  X(p,q)  &_h(p,q)  (2.22) 

Theorem  2.2  is  a  direct  corollary  of  Theorem  2.1.  This  follows 
from  the  observations  that  f(p,q)  f(-p,q)  defines  a  unitary  transforma¬ 
tion  in  L2(R3)  while  f  ■+  f*  defines  a  unitary  transformation  in  both  3C 
and  L2 (R3 ) . 
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§3 .  Transient  Free  Waves. 

The  eigenfunction  expansion  of  Theorem  2.2  is  used  in  this 

section  to  calculate  the  asymptotic  behavior  for  t  -*•  00  of  the  free 

component  u^(t,  •)  »  u(t,*).  The  principal  result  is  that  in  each  of 

the  half-spaces  R+  and  r\  where  R2  ■  {(x,y):  ±y  >  0},  uf(t,*)  is 

asymptotically  equal  to  a  wave  function  for  a  homogeneous  fluid  with 

propagation  speed  c(°°)  and  c(-°°) ,  respectively.  It  is  this  behavior 

that  motivates  the  term  "free  component"  for  u^Ct,*). 

It  will  be  assumed  that  the  total  acoustic  potential  u  satisfies 

u(t,  •)  =  Re(v(t,  •)  }  where  v(t,*)  =  exp  {-itA1^2}  h  and  h€  D(A1/,Z  )  (see 

§1).  The  corresponding  partial  waves  uk(t,*)  ■  P^  u(t,*)  with  k  >  1 

satisfy  uk(t,‘)  *  Re  (exp  (-itA1^2  )  P^h} .  This  follows  from  relations 

(2.21)  which  imply  that  Pk(h  )  *  (P^h)  .  It  follows  by  addition  that 

u  (t,*)  =  Re  {exp  (-itA1^2  )P  h}  =  Re  {v  (t,*)}  and  hence 
£  £  £ 

uf(t,#)  -  Re  (vf(t,*)}  (3.1) 

where  v^(c,#)  =  exp  (-itA1/2  )P^h  «  exp  (-itA1^2  )h^.  The  starting  point 
for  calculating  the  asymptotic  behavior  of  Uj(t,*)  will  be  the  eigen¬ 
function  expansions  of  §2.  Theorems  2.1  and  2.2  imply  the  representations 

vf(t,x,y)  *  4>+(x,y,p,q)  exp  {-i  t  A1/2  (p,q) }  h+(p,q)  dpdq  (3.2) 

f  JR3  - 

convergent  in  3C.  The  Jf-lim  notation  will  be  suppressed. 

Equation  (3.2)  gives  two  representations  of  v^  corresponding  to 
the  two  families  4>+  and  <p_.  The  calculations  below  are  based  on  the 
^-representation  which  has  been  found  to  yield  the  simplest  form  of  the 


i'tOI  11. lJp£.u 
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asymptotic  wave  function.  It  will  be  convenient  to  introduce  the  char¬ 
acteristic  functions  x+.  Xo »  and  x_  of  the  cones  C+,  C0  and  C_  in 

/v 

(p,q) -space  and  to  decompose  h_  as 


h_(P»q)  *  a(p,q)  +  b(p,q)  +  c(p,q) 


(3.3) 


where  a  »  X+  h_,  h  ■  Xo  h_  and  c  =  X_  h_.  The  corresponding  decomposi¬ 
tion  of  v^  is 


where 


va  =  exp  (-itA1//2  )ft_a 

v^  »  exp  (-itA1^2  )fi*b  > 

v£  =*  exp  (-itA1^2  )(2*c 


(3.4) 


(3.5) 


The  behavior  for  t  -*■  00  of  these  three  functions  will  be  analyzed 
separately. 

Behavior  of  va .  The  partial  wave  vfl  has  the  representation 


va(t,x,y) 


<J>_(x,y,p,q) 


exp 


(-itw+(p,q))  a(p, q)  dpdq 


(3.6) 


where 


Ui±(p,q)  -  c(±«)/|p|2  +  q2 

To  discover  the  behavior  of  v  (t,x,y)  for  (x,y)  e  R*  t  -►  »  the 

a  T 

representation  (2.6)  for  4>_<y,p,q)  -  (2q) 1/2  c(°°)  t^+(y,p,A)*  on  C+  is 
substituted  into  (3.6).  The  result  is 


1 
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va(t,x,y)  «  (2t$/T  j  exp  {i(x  •  p+yq-  taj+(p,q))}  I+(y,p,  X)"  a(p,q)  dpdq 

C+ 

C  C  °°)  I  it 

+  ( 2tt)  VI  j  exp  {i(x  •  p-  yq-  tw+(p,q)  )  }  R+(y,p,A  )  a(p,q)  dpdq 


It  is  natural  to  expect  that  in  R?  the  partial  wave  v  (t,x,y)  will 
propagate  as  t  +  ®  into  regions  where  y  is  large  and  hence  I+(y,p,X)  and 
R+(y.P>*)  ere  near  their  limiting  values.  Thus  one  is  led  to  conjecture 
that 


va(t,*)  ~  v°(t,0  +v*(t,«)  in  L2  (R^) ,  t 


(3.7) 


where 


v"(t,x,y) 


c( 


exp  {i(x  •  p+yq-  tu>+(p,q))}  a(p,q)  dpdq 


and 


1  C  ( 00  ^  I  r  it 

va(t,x,y)  °  T^-3lr  j  exp  ii(x  •  p  -  yq  -  tw+(p,q) }  R+(p,A)  a(p,q)  dpdq 


-c£3...  f 

(277)3/2  J 


exp  {i(x  •  p+yq-  tu>+(p,q))}  R+(p,A)  a(p,-q)  dpdq 


where  -C+  -  {(p,q):  (p,-q)  S  C+}  ■  {(p,q):  q  <  -a|p|).  This  conjecture 

is  proved  below.  Note  that 


v®(t,*)  -  exp  (-it  c (°°)  Ai/Z  )  h 


v^(t,*)  »  exp  (-it  c(“)  AJ/2  )  h^ 


(3.8) 


where  A0  is  the  selfadjoint  realization  in  L2(R3)  of  -A  -  -(Dz  +  Df  +  D2) 


20 


and  h  and  h1  are  the  functions  in  L2(R3)  whose  Fourier  transforms  are 
a  a  * 

A  A  ' 

h  (p,q)  “  c(°°)  a(p,q)  -  c(“)  X.  (p,q)  h_(p,q) 

>(3.9) 

h^(p,q)  -  c(°°)R+(p  A  )*a(p,-q)  -  c(“)R+(p,  A)*(l  -  x^P.q))  h_(p,-q) 

Both  are  in  L2(R3)  because  h_  6  L2(R3)  and  R+(p,A)  is  bounded,  by  (2.12). 

A  1 

Moreover,  supp  hfl  C  -  C+  and  hence  the  theory  of  asymptotic  wave  functions 
for  d'Alembert's  equation  [11,  Ch.  2]  implies  that  vx(t,  •)  ~  0  in  L,(R?) 
when  t  -*■  °°.  Combining  this  with  (3.7)  gives 

va(t,*)  -  v®(t,  •)  in  L2(R_^) ,  t  -►  °°  (3.10) 

Now  consider  the  behavior  of  va(t,x,y)  for  (x,y)  e  R3,  t  -*■  °°. 
Substituting  the  representation  (2.2)  for  $_(y,p,A)  =*  (2q) 1/2  c(<»)i[i+(y,p,A)* 
into  (3.6)  gives 

va(t,x,y)  3  exp  {i(x  •  p+  yq_  -  tco+(p,q))}  T+(y,p,A)*  a(p,q)  dpdq 

C+ 

which  suggests  the  conjecture  that 

va(t,*)  ~  v2(t,*)  in  L2(Rf),  t  -*■  “  (3.11) 

where 

va(t,x,y)  3  (2^)”^  exp  ^(x  *  P+yq_"  toj+(p,q))  T+(p,A)*  a(p,q)  dpdq 

C+ 

Now  the  mapping  (p,q)  ■+  (p,q')  -  X'(p,q)  -  (p,q_ (p,A(p,q)))  has  range 
X'(C+)  ■  R^,  Jacobian  3(p,q)/3(p,q’)  “  c2  (-®)q'/c2  (°°)q  and  satisfies 
w+(p.q)  “<*)_(p,q')-  Thus 
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v2(t,x,y) 


|  3  exP  ti(x  *  p+  yq'  -  toj_(p,q'))} 
JR+ 

X  T+(p,wf (p,q'))*  a(p,q) (c2(-«)q'/c2(“)q) 


x 


dpdq ' 


which  may  be  written 

v2  (t,  •)  -  exp  (-it  c  (— °°)  Aj/2  )  h2 
where  h2  G  L2(R3)  has  Fourier  transform 

h3(p,q')  *  c(«)  T+(p,to2(p,q))*  a(p,q(p,q'))  c2 (-“>) q '/cz (“)q 
Since  supp  h2  c  R-l  the  results  of  [11,  Ch.  2]  imply  v2(t,*)  ~  0  in 

3  »  a 

L2(R3)  when  t  -*■  00 .  Combining  this  with  (3.11)  gives 

va(t,-)  -  0  in  L2(R3),  t  -  co  (3.12) 

Analogous  conjectures  concerning  v^(t,*)  and  vc(t,*)  will  now  be 

formulated.  Only  the  main  steps  of  the  calculations  will  be  given  since 

the  method  is  the  same  as  for  v  (t,*). 

SL 

Behavior  of  vk.  Combining  (3.5)  and  (2.7)  gives 


c  (®)  r  a 

^t.x.y)  -  (  2tt)  3/2  J  exp  {i(x  •  p+yq-  tw+(p,q))}  I0(y,p,A)  b(p,q)  dpdq 


(27T)  3/2  J 


exp  {i(x  •  p- yq-  tm+(p,q))}R0(y,p,A)  b(p,q)  dpdq 


which  suggests  the  conjecture 


vfc(t,‘)  ~  v°(t,*)  +v£(t,')  in  L2(R2),  t 


-►  00 


(3.13) 
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where 


v°(t,*)  =  exp  (-it  c(°°)  Aj/2  )  h. 


(3.14) 


v£(t,*)  =  exp  (-it  c(«)  Aj/2  ) 
and  and  are  the  functions  in  L2(R3)  whose  Fourier  transforms  are 
^(p.q)  «  c(")  b(p,q)  =  c(°°)  XoCP*^)  h_(p,q) 


h£(p,q)  =  c  (°°)  R0  (p,A)*  b(p,-q)  =>  c(°°)  R„(p,X)*  (l-x0(P»q))  h_(p, 


t>  (3.15) 

q) 


Since  supp  h^  C  -  CQ  C  R3,  v^(t,*)  ~  0  in  L3(R_^)  when  t  by  [11, 

Ch.  2].  Combining  this  with  (3.13)  gives 


vb(t,«)  ~  v®(t,-)  in  L2(r|),  t  -*• 


Similarly  (3.5)  and  (2.3)  imply  the  representation 


(3.16) 


vb(t,x,y) 


(2  IT) 
and  since 


exp  (i(x  •  p  -  taj+(p,q))  }  T0(y,p,X)"  exp  (gy/H)  b(p,q)  dpdq 
Co 


11m  T0(y,p,A)  exp  (gy/H) 

y  ♦■—oo 


one  expects  that 


vfe(t,  •)  '  0  in  L2  (R3) ,  t  -*• 


(3.17) 


Behavior  of  vr .  Combining  (3.5)  and  (2.2)  gives 
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vc(t,x,y) 

exp  { i(x  •  p-  yq^.-  tto_(p,q))}  T_(y,p,A)*  c(p,q)  dpdq 
which  suggests  that 

vc(t,*)  ~  v*(t,*)  in  L2  (R*) ,  t  ->•  «  (3.18) 

where 


v^(t ,x,y) 


c(-°°) 
(2tt)  3/2 


exp  { i(x  •  p  -  yq+  -  t(q_(p  ,q) )  }  T_(p , A)  c(p,q)  dpdq 


The  mapping  (p,q)  (p,q')  “  X"(p,q)  *  (p.-q+Cp.A^.q)))  has  range 

X"(C_)  ■*  -  C+,  Jacobian  3(p,q)/3(p,q')  =  c2  (°°)qVc2  (-®)q  and  satisfies 
w_(p,q)  *  w+(p,q').  Thus 

v*(t,«)  -  exp  (-it  c  (°°)  Aj/2  )  h^ 

where  h1  has  Fourier  transform 
c 

h*(p,q')  *  c(-o°)  T_(p,oi2(p,q’))*  c(p,q(p,q’))  c2  (“)  q '/c2  (-«)q 

Moreover,  supp  h*  c  -  C+  C  R2  and  hence  one  expects  v*(t,*)  -  0  in 
l2(r+)»  t  Combining  this  with  (3.18)  gives 


vc(t,  •)  ~  0  in  Lz  (R+) ,  t 


(3.19) 


Finally,  combining  (3.5)  and  (2.2)  gives 


(t,x,y)  -  (2tt)^  exp  {i(x  •  p+ yq  -  tw_(p,q))  }  I_(y,p,A)  c(p,q)dpdq 
C 


(2ir) 


$r 

f 


exp  (i(x  •  p-  yq-  tu>_(p ,  q)  )  }  R_(y ,  p ,  A)  c(p,q)  dpdq 


which  suggests  that 


vc(t,*)  ~  v®(t,')  in  L2  (R*) ,  t  (3.23) 

The  asymptotic  behavior  of  v^(t,*)  for  t  »  is  obtained  from 
the  three  cases  analyzed  above  by  superposition,  equation  (3.4).  Thus, 
equations  (3.10),  (3.12),  (3.16),  (3.17),  (3.19)  and  (3.23)  imply 

+  vb(t’,:)  in  MR+> 

Vf  (t,  •)  -  <  >  t  -*•  « 

i  v“(t,*)  io  Lj(Ri)  , 

Combining  this  with  the  definition  of  v®,  v£  and  v®,  equations  (3.8), 
(3.9),  (3.14),  (3.15),  (3.21)  and  (3.22)  suggests 

Theorem  3.1.  For  all  h  €  TC  let  v|(t,*)  be  defined  by 
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v?(t,x,y)  = 


where  h  and  h 
given  by 


exp  (-it  c(°°)  A*/2  )  h+(x,y),  (x,y)  e  R+ 


(3.24) 


{  exp  (-it  c ( — °°)  A^2  )  h~(x,y)  ,  (x,y)  €  R3  J 

are  the  functions  in  L2(R3)  whose  Fourier  transforms  are 


h+(p,q) 


c(°°)  h_(p,q) ,  (p  ,q)  £  R^ 
0  ,  (p , q)  €  R3 


' 

> 


and 


h  (p,q) 


0  ,  (p,q)  e  R+ 

c(-°°)  h_(p,q),  (p , q)  e  R3 


> 


Then 


lim  II  v  ( t ,  *  ) 

£-K»  1 


Vj(t,  * ) II =  0 


(3.25) 


(3.26) 


(3.27) 


Proof  of  Theorem  3.1.  The  decomposition  (3.3)  will  be  used  for 
the  proof.  Moreover,  for  brevity,  only  the  asymptotic  equality  (3.7) 
for  v  (t,«)  will  be  proved.  The  five  remaining  cases,  namely  (3.11), 
(3.13),  (3.17),  (3.18)  and  (3.20),  can  be  proved  by  the  method  used 


for  v  . 
a 

As  a  first  step,  (3.7)  will  be  proved  for  the  special  case  where 
a(p>q)  €  c”(C+) .  The  general  case  will  then  be  proved  by  using  the  fact 
that  CQ(C+)  is  dense  in  L2 (C+) . 

When  a(p,q)  6  C  (C.)  the  integrals  defining  v  ,  v°  and  v1 

u  "T"  3  3  3 

converge  point-wise,  as  well  as  in  7f,  ana  one  has 
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v  (t,x,y)  -  v®  ( t ,x, y)  -  v* (t ,x, y) 


c(”) 

(2TT)  3/2 


exp  (ix  •  p) w(t,y ,p)  dp 


(3.28) 


where 


w(t,y,p)  =  w°(c,y,p)  +  w1 (t,y,p) 


(3.29) 


w°(c,y,p)  = 


a  P 


exp  {i(yq  -  tco+(p,q)  )  }(I+(y,p,A)  -  1)  a(p,q)dq  (3.30) 


and 


w^t.y.p)  =*  I  exp  {i(-yq  -  tw+(p,q)) } (R+(y,p,A)  -  R+(p,A))  a(p,q)dq 

a|  p  | 

(3.31) 

ralp|  * 

exp  { i(yq  -  tw+(p,q)) }(R+(y,p,A)  -  R+(p,X))  a(p,-q)  dq 


Parseval's  formula  In  L2(R2),  applied  to  (3.28),  gives 


2  |va(t,x,y)  -  v°(t,x,y)  -  v*(t,x,y) | z  dx  =  — ^  |w(t,y ,p) | 2  dp 


Integrating  this  over  y  >  0  gives 


,va(t,-)  -  v«(t.*)  -  =  j^Tr  l»(t.*)lIi2(R,) 


(3.32) 


c(”) 


-  (2ir) 


T7T 


,"w0(t.*>llL2(R3)  +  Bwl(t’*)'lL2(R|) 


by  (3.29)  and  the  triangle  inequality. 

The  estimate  (3.32)  implies  that  to  prove  (3.7)  it  is  sufficient 
to  prove  that  w°(t,*)  ■*  0  and  w1(t,*)  -*■  0  in  L2 (R^)  when  t  -*•  °°.  To  this 


end  the  integrals  in  (3.30)  and  (3.31)  will  be  estimated  by  the  method 
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of  stationary  phase  as  formulated  in  the  Appendix.  To  apply  the  method 
define 

r  =  / t2+y2 ,  t  =  r  sin  ot,  y  =  r  cos  a  where  0  <  a  <  tt /  2 


and 


9(p,q,a)  =  q  cos  a  -  ui+(p,q)  sin  a  (3.33) 

Then  the  integrals  in  (3.30)  and  (3.31)  have  the  form  (A.l)  with 
s  =  q  (m  =  1)  and  £  =  (p,a)  (n  =  3).  Moreover, 


where 


0(p,q,a) 


cos  a  -  U(p,q)  sin  a 


U(p,q)  -  c(“)q//|p |2+q2 

is  the  group  speed  associated  with  the  dispersion  relation  w  =  a>+(p,q). 
Clearly  0(p,q,a)  ■  0  if  and  only  if 

y/t  ■  ctn  a  *  U(p,q) 

Note  that  ctn  a  >  0  for  0  <  a  <  tt/2.  In  the  case  of  (3.31),  U(p,q)  <  0 
on  the  interval  of  integration  and  there  are  no  points  of  stationary 
phase.  In  the  case  of  (3.30),  a|p|  <  q  <  00  and  hence  0  <  U(p,q)  <  c(°°). 
Thus  there  is  a  unique  point  of  stationary  phase  if  0  <  ctn  a  <  c(°°),  or 
a0  <  a  <  tt/2  where  ctn  ag  *  c(«),  and  no  point  of  stationary  phase  if 
0  <  a  <  a0 .  The  point  of  stationary  phase  is 

q  =*  x(p,oi)  *  |  p  |  ctn  a/(c2(°»)  -  ctn2  a) 1/2  (3.34) 


when  0  <  ctn  a  <  c(°°).  The  corresponding  approximation  to  w°(t,y,p)  is, 


from  (A. 11)  , 


w  (t,y,p) 


(3.35) 


X(t,y)  (2ir)1/z  g(r,p,q,a)  exp  {i(yq-  tw+(p,q))  -  tt/4>/ ( t  U  (p,q))1/z 


where  x(t,y)  Is  the  characteristic  function  of  the  set 
{ (t,y) :  0  <  y  <  c(°°)t}. 


g(r,p,q,a)  =  (I+(r  cos  a,p,A(p,q))  -  1)  a(p,q)  (3.36) 

Uq(p,q)  =  c(“)|p|2/(|p|2+q2)3/2 

and  q  is  given  by  (3.34)  with  ctn  a  =  y/t.  These  results  and  Theorems 
A. 1  and  A. 2  imply  the 

CO  * 

Lemma.  For  every  a(p,q)  €  C0(C+)  there  exists  a  constant 
M  =  M(a)  such  that 


|w°(t,y,p)  -  w°°(t,y,p)  |  <  M/(t2  +  y2)3/l* 


(3.37) 


for  all  t  >  0,  y  >  0  and  p  G  R2  .  Similarly,  for  each  positive  integer  k 
there  exists  a  constant  =  Mk(a)  suc^  that 

[w1 (t,y,p) |  <  Mk/(t2  +  y2)k/2  (3.38) 

for  all  t  >  0,  y  >  0  and  p  G  R2 . 

Proof  of  the  Lemma.  The  hypothesis  a  e  C*(C+)  implies  that 
supp  a  ■  K0  is  a  compact  subset  of  C+.  Let  K  and  Kj  denote  its 
orthogonal  projections  on  the  q-axis  and  the  p-plane,  respectively.  Then 
the  functions  w  (t,y,p),  w  (t,y,p)  and  w  (t,y,p)  are  all  zero  for 
p  €  R2  -  Kj .  Hence  it  is  enough  to  verify  (3.37)  and  (3.38)  for  p  €  Kj . 
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ft 


f 


ft 


i 


* 


ft 


ft 


ft 


» 
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consider  first  the  integral  in  (3.31).  It  has  the  form  (A.l) 
with  q  6  0  =  R_  (the  negative  real  numbers),  (p,a)  e  O' 

=  (R2  -  { 0 } )  x  (0,tt/2)  and  g  given  by 

g(r,p,q,a)  *  [R+(r  cos  a,p,A(p,q))  -  R+(p, A(p ,q) ) ] *a(p ,-q)  (3.39) 

Hypothesis  (A. 2)  is  evident  from  (3.33),  while  (A. 3)  follows  from  (3.39) 
and  well-known  properties  of  hypergeometric  functions.  Moreover, 
supp  g(r,*,*,a)  C  K  for  all  r  >  0  and  0  <  a  <  tt/2  and  hence  (A. 4)  holds. 
The  uniform  estimates  (A. 6)  may  be  derived  from  the  integral  representa¬ 
tions  of  the  hypergeometric  functions  -  see  [3,  pp.  77-79].  Finally, 
0(p,q,a)  has  no  points  of  stationary  phase  for  g  £  K  C  R_.  Thus  Theorem 
A.l  is  applicable  to  w^t.yjp)  and  (3.38)  follows. 

Next  consider  the  integral  in  (3.30).  For  0  <  a  *  ctn-1(y/t) 

<  a0  there  are  no  points  of  stationary  phase  q  £  K.  Moreover, 

00 

w  (t,y,p)  =  0  on  this  set.  Hence  Theorem  A.l  is  applicable  and  implies 
that  (3.37)  holds  for  y  >  c(°°)  t.  On  the  other  land,  for  aQ  <  a  <  tt/2 
or  0  <  y  <  c(°°)t  there  is  exactly  one  point  of  stationary  phase  q  given 
by  (3.34)  and  hence  (A. 8)  holds.  Moreover, 

6q(p,q»a)  *  -Uq(p,q)  sin  a  ¥  0  for  q  £  K,  (p,a)  £  O' 

which  implies  (A. 9).  Finally,  (A. 10)  follows  from  properties  of  the 
hype rgeome trie  functions,  as  in  the  preceding  cases.  Thus  Theorem  A. 2 
is  applicable  and  implies  that  (3.37)  holds  for  0  <  y  <  c(«>)t. 

Proof  of  Theorem  3.1  (concluded).  The  lemma  implies  that  if 
a  €  C0  (C+)  then  w^t,*)  and  w1(t,*)  tend  to  zero  in  L2  (R^,)  when  t  ■*  00 
and  thereby  proves  Theorem  3.1  in  this  case.  Indeed,  integrating  (3.38) 
gives 


ft 
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lwl(t,')llL2(R^)  =  J 


<  |kJ 


|w‘ (t,y,p) | 2  dpdy 


(t2  +  y2)~K  dy  »  t 


(3.40) 


where  | Kj  j  is  the  measure  of  Kj  and  *  c£  | K x  |  /  (1  +  £2)  ^  d£.  Now 
consider  w°.  The  lemma  implies  that 


w°(t,y,p)  =>  w  ( t ,  y ,  p)  +  q(t,y,p) 


where 


|q(t,y,p)  |  <  M/(t2  +  y2)3/“ 


Thus 


« w°  (t,  *)I1L^  (R^,>  <  II  w  <c  » *  >  II  Lz  (r_3)  +  flq(t,  _ )  Kl2  (R^) 


(3.41) 


Moreover, 


'■>(t-->'L2(E|)  2C,/t’ 

by  the  calculation  used  to  prove  (3.40)  with  k  *  3/2.  It  remains  to 
show  that  #w°°(t,  ')HL  0  when  t  -*•  °°.  Direct  calculation  using 
(3.35)  and  (3.36)  gives 


]w  (t,*)l 


2tt 


L2 (R^) 


c(“»)t 


c(“)t 


[  |w”(t,y. 


p)|2  dpdy 


(3.42) 


l+(y.p»*(p.q))  -  l|2  i a(p  ,q) I  2  t_1  U  (p,q)_1  dpdy 


where  q  *  x (p , ctn-1 (y/t) )  is  given  by  (3.34).  Since  (3.34)  defines  the 
solution  q  of  y/ t  =  U(p,q)  it  is  natural  to  make  the  change  of  variable 
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y  =  t  U(p,q),  dy  =  t  U^(p,q)  dq  (p  fixed)  in  (3.42), 


The  result  is 


lw  (t,,)llL2(R|) 


=  2tt 


(3.43) 


[  !l+(t  U(p,q) ,p,A(p,q))  -  lj2  |g(p,q)j2  dqdp 

Kj  J  K 


since  supp  g  C  Kj  x  K.  Finally  the  representation  I+(y,p,X) 
=  (2/H)1/2exp  (iq+y)  02(y,p,X)  implies  that 

I+(t  LT(p,q)  ,p,X(p,q))  -*  1  v/hen  t  ■>  00 


uniformly  for  (p,q)  t:  tC1  x  K.  It  follows  by  Lebesgue’s  dominated 
convergence  theorem  that  H  w°°(t,  *)flT  ,_3n  =  0(1),  t  In  fact,  a 

careful  estimate  of  the  hypergeometric  function  <)>u  gives  the  stronger 
estimate 


I!  w°°( t ,  • ) 11  Lz  (R3  -j  =  0(e  Mt)  ,  t  -*•  °° 

where  p  ■  p(a)  is  a  positive  constant. 

The  arguments  given  above,  applied  to  v  ,  v,  and  v  ,  show  that 

3D  C 

the  conclusion  (3.27)  of  Theorem  3.1  holds  for  all  h  such  that 
h_  e  C0(C+  U  C0  U  C_) .  Moreover,  this  set  is  dense  in  L2(R3)  and  hence 
ft_C0 (C+  u  C0  U  C_)  is  dense  in  m  ?^K  by  Theorem  2.2.  These  facts  can 
be  used  to  extend  (3.27)  to  all  h  S  JC,  provided  the  mappings 
U(t) :  Xf  *  L2(R3)  and  U0(t):  3ff  -*■  L2(R3)  defined  by  U(t)hf 
*  exp  (-it  A1^2  )hj  and  U0(t)h^  »  v£(t,*)  are  uniformly  bounded  for  all 
t  e  R.  This  density  argument  has  been  given  in  many  places;  see,  for 
example,  [11,  Ch.  2]  or  [13,  p.  260].  To  prove  the  uniform  boundedness 
note  that  the  L2(R3)  and  X  norms  are  equivalent: 
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c(y0)  Jhljf  <  9h«L  <  c(-°°)  a  till  ^ 

and  hence 

IU(t)hflL  (r3)  <  c(-»)  !U(t)hf  11^  =  c(-»)  Ihjlj,. 
for  all  t  6  R.  Similarly,  (3.24),  (3.25)  and  (3.26)  imply  Che  estimates 


“Uo(t)hfSL2(R3) 


»v?(t, ■ 


)t 


M*3) 


-  (,Vf(t’,),L2(Rj)  +  “v?(t**)"L2(R2))1/2 


^(“hlL2(R3)  +  “h"»i2(R3))1/2 


-  ^+“L(R^)+“^L(R3)>1/2 


-  (c2(»)  .hjl*^  +  C>(-)  lfiji>(Rl))l,i 

<  c(-«)  «h.lL2(R3)  =  c(-oo)  lhflK 

for  all  t  €  R.  This  completes  the  proof  of  Theorem  3.1. 

Theorem  3.1  implies  corresponding  asymptotic  estimates  for  the 
free  component  uf(t,*)  *  Pf  u(t,*)  *  Re  {vf(t,*)}  of  the  acoustic 
potential  u(t,’).  Indeed,  if  u£(t,*)  is  defined  by 

u°(t,x,y)  -  Re  {v°(t,x,y)>  (3.44) 


O 


then  Theorem  3.1  and  the  elementary  inequality  |Re  z|  <  |z|  imply 
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Corollary  3.2.  For  all  h  G  3C  one  has 


lim  ||uf(t,*)  -  *  0 

t-"» 


(3.45) 


If  the  Initial  state  h  has  derivatives  in  3f  then  u^(t,*)  and 
u°(t,*)  have  the  same  derivatives  in  K  and  (3.45)  can  be  strengthened  to 
include  these  derivatives.  In  particular,  one  has 

Corollary  3.3.  For  all  h  e  l£(R3)  =  D(A1/2  )  one  has 


lim  JD,  u.(t, • )  -  D.  u^(t, *)ljf  =  0,  j  =  0,1, 2, 3 


(3.46) 


t-KO 


where  D0  =  3/3t,  Di  =  3/3xlf  D2 =  3/3x2  and  D3  *  3/3y.  Note  that  (3.46) 
is  equivalent  to  convergence  in  energy: 


lim  E(uf  -  u®,R3,t)  =  0 
t-**> 

Corollary  3.3  can  be  proved  by  applying  the  method  of  this 
section  to  the  derivatives  vf(t,*)  (j  *  0,1, 2, 3)  which  are  given  by 

integrals  of  the  same  form  as  (3.2).  Detailed  proofs  for  the  case  of 
the  Pekeris  profile  were  given  in  [13]. 
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§4.  Transient  Guided  Waves. 

The  asymptotic  behavior  for  t  -*■  “  of  the  guided  component 

u  (t,  •)  *  P  u(t, •)  is  derived  in  this  section,  u  (t,*)  is  a  sum  in  Jf 
S  o  S 

of  mutually  orthogonal  partial  waves  u^(t,*)  *  u(t,*) 

*  Re  (vk(t,*)},  k  >  1.  The  starting  point  for  the  analysis  is  the 
integral  representation 


vk(t,x,y) 

where 


j  ^k(x,y,p)  exp  {-it  o^Cjpl)}  hk(p)  dp 
Qk 


(4.1) 


\(p)  =  J  3  \(x,y,p)*  h(x,y)  c~Z (y)  dxdy 

and  the  integrals  converge  in  Jf  and  Lz  (&k) ,  respectively.  The  integral 
in  (4.1)  can  be  written 


vk(t,x>y)  -  ^ 


“k 


exp  (i(x  •  p  -  tc^dpl))}  4<k(y»p)  hk(p)  dp  (4.2) 


This  is  an  oscillatory  integral  that  can  be  estimated  by  the  method  of 
stationary  phase  of  the  Appendix  when  hk  €  (^(f^) .  To  apply  the  method 


define 


r  =*  /tz  +  j  x | 2 

t  1  1  (0,  x,  »  r  x2  -  r  C2 
5  -  (5o.Cx.5a)  6  S2  c  R3 

where  S2  denote  the  unit  sphere  in  R3 .  Then  (4.2)  takes  the  form  A.l 


bi^vnK-wOT  HusEli 


'■  4  - 


36 


with  m  =  2,  n  *  3 


0k(p.£.y)  =  SiPi  +  C2p2  -  €oVIpI>  (4'3) 

gk(r.P»C.y)  *  (2tt)'1  ^k(y,p)  hk(p)  (4.4) 

Note  that  8^  is  independent  of  y  and  gk  is  independent  of  r  and  £•  They 
have  been  written  in  the  forms  (4.3)  and  (4.4)  to  emphasize  the  applica¬ 
bility  of  the  results  of  the  Appendix. 

The  phase  function  (4.3)  has  a  point  of  stationary  phase  if  and 

only  if 


Uk(jpi)p/|p|  =  x/t  (4.5) 

where 

uk(|p|)  3  D|p|  V|p|) 

is  the  group  speed  associated  with  the  dispersion  relation  w  *  cok ( | p | ) . 
The  analysis  will  be  simplified  by  the  following 

Proposition.  If  the  parameters  defining  c(y),  equation  (1.1), 

satisfy 


8K  >  M  (4.6) 

then  each  of  the  functions  Uk(|p|),  k  -  1,2,***,  is  a  strictly  decreasing 
function  on  pk  <  |p|  <  ®  such  that  Uk(pk)  »  c(®)  and  Uk(|p|)  c(y0)  when 

Ip!  -  *. 

A  proof  of  the  proposition  may  be  derived  from  the  parametric 


9 


O 


4 


representation  of  the  dispersion  relation  [4].  The  calculation  is  too 
long  to  give  here. 
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For  simplicity  of  presentation  condition  (4.6)  is  assumed  in  the 
remainder  of  this  paper.  If  8K  <  M  then  U£(|p|)  may  have  a  finite 
number  of  zeros.  This  complicates  the  form  of  the  asymptotic  wave  func¬ 
tions  but  is  otherwise  tractable.  In  the  case  of  the  Pekeris  profile, 
treated  in  [13],  the  functions  U^(|p|)  have  exactly  one  zero  on 
I P I  >  Pk- 

By  calculating  the  Hessian  0^  one  can  show  that 
r2  |det  6£(p,C,y)|  -  t2  Ufc( | p | )  |U£( | p | ) | / | p | 

and 

sgn  0£(p,e.y)  -  0 

In  particular,  each  point  of  stationary  phase  is  non-degenerate,  by  the 
Proposition.  Moreover,  (4.5)  implies  that  9^  has  a  point  p  of  stationary 
phase  if  and  only  if 

|x|/t  -  Uk( | p  j )  (4.7) 

Now  (4.7)  has  a  unique  solution  |p|  if  |xj/t  lies  in  the  range  of 
Uk(|p|);  that  is  if 

c(y0)  <  |x|/t  <  c(®)  (4.8) 

The  point  of  stationary  phase  is  then  given  by 

p  -  Q^( | x j /t)x/ | xj  for  x  €  R2  -  {0},  t  >  0  (4.9) 

where  Qk  is  the  inverse  function  to  U^.  It  is  a  monotone  decreasing 
function  that  maps  (c(y0)  ,c(°°))  onto  (Pk>°“)  •  By  Theorem  A. 2  the  point 
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(4.9)  makes  a  contribution 

„  |pi1/2exp  {  i(  jx  |  |  p  |  -  tui.  (  |  p  |  ) )  )  \(y,p)  h,(p) 

v  (t,x,y,p)  »  -  (4.10) 

k  t  {uk(|Pj)  lu'dpDl}^ 

to  the  integral  in  (4.2)  when  (t,x)  satisfies  (4.8).  For  |x|/t  outside 
the  interval  (4.8)  there  is  no  point  of  stationary  phase.  Thus  the 
stationary  phase  approximation  to  vk(t,x,y)  is  given  by 

vk(t,x,y)  =  x(|x|/t)  v”(t,x,y,Qk( |x|/t)x/ |x|)  (4.11) 

where  x  is  the  characteristic  function  of  (c (y0 )  ,c(°°) )  and  one  has 

Theorem  4.1.  For  all  h  €  K  such  that  hk  6  CQ(f2k)  there  exists  a 
constant  C  =  Ck(h)  such  that 

|vk(t,x,y)  -  v"(t,x,y)|  <  C/t2  (4.12) 

for  all  t  >  0,  x  €  R2  -  {0}  and  y  6  R. 

This  result  can  be  proved  by  application  of  Theorems  A.l  and  A. 2. 
The  proof  is  similar  to  that  of  the  Lemma  in  §3  and  will  not  be  recorded 

here.  If  hk  is  not  a  smooth  function  then  Theorems  A.l  and  A. 2  are  not 

applicable  and  the  estimate  (4.12)  may  fail.  However,  the  definition 

(4.10) ,  (4.11)  is  applicable  to  all  h  S  Jf.  More  precisely,  one  can  prove 

Theorem  4.2.  For  all  h  e  3C,  all  t  >  0  and  k  »  1,2, •••  one  has 

v*(t,‘)  6  Jf, 

and 

avk(c’  * ) * 7f  *  “  ^k^Jf 
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Moreover,  Che  mapping  t  -*■  v^(c,*)  is  continuous  from  R+  to  JC  and 

lim  llvk(t,*)  -  =  0  (A.  13) 

t-«0  *" 

The  proofs  of  these  properties  are  the  same  as  those  for  the 
Pekeris  profile,  given  in  [13],  and  are  not  reproduced  here.  On  defining 


u"(c,x,y)  -  Re  {vk(t,x,y)} 


one  has 

Corollary  4.3.  For  all  h  <=  K  and  k  =  1,2,3,*»* 


lim  Juk(t,*)  -  uk(t, •)#,/•  *  0 


C-*®0 


X 


If  h  G  Lj(R3)  then  uk(t,*)  6  L^R3)  and  asymptotic  wave  functions 
for  the  first  derivatives  of  uk  can  be  constructed.  Indeed,  if 

~  OO 

hk  S  co then  the  first  derivatives  of  vfc  are  given  by 


Dt  vk(t,x,y) 
Dj  vk(t’x’y) 
Dy  vk(t’x’y) 


2tt 


% 


exp  {i(x  •  p  -  tou^ ( |  p | ))  )  <l'k(y,p)(-i<i>k(|p|))  hk(P>  dP 


2tt 


exp  { 1  (x  •  p  -  ta^ ( |  p |  ) )  }  ^k(y,p)(iPj)  hk(p)  dp  (j  -  1,2) 


_i_ 

2tt 


^k 


exp  {i(x  •  p  -  ta^dpl  ))}  D  (]<k(y,p)  \(p)  dp 


These  integrals  have  the  same  form  as  the  integral  (A. 2)  for  vk(t,x,y). 
The  corresponding  asymptotic  wave  functions  are  defined  by 


vko (t,x,y  ,p)  =*  < - iu>k ( |  p | ) )  v"(t,x,y,p) 
v"j(t,x,y,p)  -  (ipj)  v“(t,x,y,p)  (j  -  1,2) 
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! 


b 


v”3(t,x,y,p)  =  Dy  v“(c,x,y,p) 

and 

v“j(t,x,y)  =  X(lx)/t)  v“j (t,x,y,Qk(Jx|/t)x/Jx|)  (j  =  0,1, 2,3) 

The  analogue  of  Theorem  4.2  is 

Theorem  4.4.  For  all  h  6  L^CR3),  all  t  >  0  and  k  =  1,2,3,*** 

one  has 


vk0(t,*)  6  Jf,  vkj(t,*)  €  L2(R3),  j  =  1,2,3 


#vko  +  ^  l,Vkj^t’’')[lL2(R3)  ~  2  IIAl/  hkllJf 


and 


lim  ID.  Vj(t ,  • )  -  vk-(t,*)U7f  =  0  (j  =  0,1, 2, 3) 
£->O0  -J  d 


Finally  if 


uk^.(t,x,y)  *  Re  {v*..  (t  ,x,y)  }  (j  =  0,1, 2,3) 
then  one  has 

Corollary  4.5.  For  all  h  e  L2(R3)  and  k  =  1,2,3, ••• 
lim  II D.  uk(t,*)  -  u"  ( t , * ) II ^  *  0 
Theorem  4.4  may  be  proved  by  the  technique  used  for  the  Pekeris 


profile  in  [13] . 


.'«is.K ...  .■» 
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§5.  The  Asymptotic  Distribution  of  Energy  for  Large  Times. 

The  total  energy  of  the  acoustic  field  u(t,x,y),  given  by  (1.16), 
is  constant  for  t  >  T.  The  same  is  true  of  the  partial  waves  u^,  u^  and 
u^,  k  =  1,2,*“.  Moreover,  the  eigenfunction  expansion  of  [4]  implies 
that  {Pj,Px ,P2 , • • • }  is  a  complete  family  of  orthogonal  projections  that 
reduces  A.  It  follows  that 


JA1/2  hll^  =  II Al/2  hf||*.  +  l  |Al/2hkl2Jf  (5.1) 

k=l 

which  may  be  interpreted  as  an  energy  partition  theorem.  The  energies  of 
the  partial  waves  can  be  calculated  from  the  initial  state  h,  or  source 
function  f(t,x,y),  by  means  of  the  eigenfunction  expansion  theorem.  The 
relationship  between  h  and  f  is  given  by  (1.9).  Hence 

h_(p,q)  =  i  A  1/2  (p ,  q)  f_(-Al/2  (p,q),p,q) 


where 


f_(w,p,q) 


exp  (-iiox)  <p  (x,y,p,q)  f(x,x,y)  c  z(y)  dydx 


and  similarly 


-1 


\(P)  -  i  <  !  P  I )  f  I  P I  )  .P) 


where 


f ,  (oj  ,  p ) 


exp  (-iu)T )  i^,(x,y,p)  f(x,x,y)  c"2(y)  dydx 
>  3  K 
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It  follows  that  the  partial  energies  are  given  by 


E(uf ,R3 ,t) 


1/2  h  «  2 


|A1/Z  h-| 


f'3C 


|f_(-X1/2  (p ,q) »p ,q) | 2  dpdq,  t  >  T 


and 


E(uk,R3,t)  = 


1Al/2  Vic 


lM-(M|p|  )  .p)  I  2  dp,  t  >  T 


The  theorems  of  §3  and  §4  make  possible  the  calculation  of 
asymptotic  energy  distributions  in  bounded  and  unbounded  subsets  of  R3 . 
Only  the  principal  results  are  formulated  here  since  the  proofs  are  the 
same  as  for  the  case  of  the  Pekeris  profile  given  in  [13]. 

The  notation 

E°°(u,K.)  ■  lim  E(u,K,t) 

f+<» 

will  be  used  whenever  the  limit  exists.  A  first  result  is  the 
transiency  of  all  waves  with  finite  energy  in  an  Epstein  duct: 

OO  . 

E  (u,K)  =  0  for  all  compact  sets  K  c  R3 

Next  let  C+  (resp.  C~)  denote  a  cone  in  R^  (resp.  R3) .  Then  Corollary 
3.3  and  the  results  of  [11]  imply 


E°0(uf,C±)  -  c2(±«)  (  +  ( |  p  |  2  +  q2)  |  h  (p,q)  |  2  dpdq 
>C~ 

"  *  |f  (-Al/2  (p,q)  ,p,q)  |  2  dpdq 

JC" 


Now  define  the  cone  C_  by 


A3 


C_  =  { (x,y) :  f  y  -  y  o |  <e|x|},£  >0 

Then,  in  contrast  to  the  preceding  result,  one  has 

e^Cu^.C^ )  *  e(uk,R3,t)  =  JA1/2  h^II^.  for  every  e  >  0 

and  for  k  *  1,2,3, ••*  (see  [13,  Theorem  5.5]). 

Finally,  if 

S  *  {(x,y):  x  e  R2  and  y5  <  y  <  y2l 
then  it  can  be  shown  by  means  of  Corollary  A. 5  that 

00  (•  f  ry2 

e  (u  ,s)  -  | f.  (-to,  (|p|) ,p) | 2  4c(y,p)  c  2(y)dy  dp 

\  lJyx  k 


I 
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§6.  Concluding  Remarks.  Examination  of  the  proofs  in  [4]  and  in  this 
paper  reveals  that  the  methods  should  apply  to  a  large  class  of  sound 
speed  profiles  c(y).  As  remarked  in  [4],  the  principal  obstacle  to  such 
a  generalization  is  the  lack  of  information  about  the  p-dependence  of 
the  eigenfunctions  and  eigenvalues.  In  the  cases  of  the  Pekeris  and 
Epstein  profiles  this  information  was  obtained  from  explicit  construc¬ 
tions  of  the  eigenfunctions  in  terms  of  known  functions. 


i-'rtLCLuUG  irr.-S.  bJu*i«K-WOT  FliJCA 
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» 


> 
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Appendix.  Estimates  of  Oscillatory  Integrals  with  Parameters. 

The  method  of  stationary  phase  provides  asymptotic  estimates  for 
r  -*■  00  of  oscillatory  integrals  of  the  form 

I(r)  *  [  exp  (ir  9(s)}  g(s)  ds,  0  c  Rm 

The  method  is  needed  in  §3  and  §4  above  to  estimate  the  oscillatory 
integrals  in  equations  (3.30),  (3.31)  and  (4.2).  In  (4.2)  the  integrand 
g  contains  a  parameter  y  in  addition  to  the  variables  of  integration.  In 
(3.30)  and  (3.31)  both  0  and  g  contain  parameters  and,  in  addition,  g 
contains  the  large  parameter  r.  In  both  cases  estimates  that  are 
uniform  in  the  parameters  are  needed.  This  suggests  the  study  of 
oscillatory  integrals  of  the  form 


I(r,F.) 


i  exp  {ir  0(s,£) }  g(r,s,£)  ds 

J0 


(A.l) 


where  r  >  0,  s  E  Rm  and  £  €  Rn.  9(s,Q  is  a  real-valued  phase  function 
and  it  is  assumed  that  there  are  open  sets  0  C  Rm  and  O'  C  Rn  such  that 

Dg  9(s,£)  €  C(0  x  O')  for  all  multi- indices  6  (A. 2) 

g(r,s,£)  €  C(R+  x  0  x  O’)  for  all  multi- indices  6  (A. 3) 

where  R+  denotes  the  positive  real  numbers.  Moreover,  it  is  assumed 
that  there  is  a  compact  set  K  C  0  such  that 

supp  g  C  R+  x  R  x  O’  (A. 4) 


irtL  '.Luj..O  bL*uiK-i*Or 
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Estimates  of  I(r,£)  for  r  -*•  00  are  sought  which  are  uniform  in  £  on  compact 
subsets  of  O'.  For  large  values  of  r  the  exponential  in  (A.l)  is  highly 
oscillatory  except  near  critical  points  of  the  phase  function  9(s,£). 

Two  estimates  of  I(r,£)  are  given  here,  corresponding  to  the  cases  of  no 
critical  points  and  one  critical  point,  respectively.  The  first  case  is 
formulated  as 

Theorem  A.l.  Assume  that 

Vg  0(s ,£)  =  (38/3s1,***,36/3sm)  0  for  all  (s,£)  e  K  x  O'  (A. 5) 

Moreover,  assume  that  for  each  compact  set  K ’  C  O' ,  each  r0  >  0  and  each 
positive  integer  k  there  exists  a  constant  M  =  M(K,K' ,r  ,k)  >  0  such  that 

£ 

|D  g(r,s,£)|  <  M  for  all  r  >  r  ,  s  e  K,  £  e  K'  and  |6|  <  k  (A. 6) 

Then  there  exists  a  constant  C  =  C(K,K' ,r0 ,k,g)  >  0  such  that 

|l(r,£)|  5  C  r_k  for  all  r  >  r0  and  E,  €  K'  (A. 7) 

In  the  second  case  considered  here  0(s,£)  has  a  unique  non¬ 
degenerate  critical  point  s  =  x(£)  for  each  £  e  O'.  It  is  formulated  as 
Theorem  A. 2.  Assume  that  there  is  a  function  T  €  C  (O' ,0)  such 
that,  for  all  £€(?', 

Vg  9(s,£)  *  0  if  and  only  if  s  *  t(£)  (A. 8) 

Moreover,  assume  that  the  Hessian  9"(s,£)  »  (320(s,£)/3s.3s,  )  satisfies 

J  K 

det  0,,(t(C),5)  i  0  for  all  Z  e  O'  (A. 9) 

In  addition,  assume  that  for  each  compact  set  K'  c  O'  and  each  r„  >  0 
there  exists  a  constant  M  »  M(K,K',r0)  such  that 
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-  i>. 


|Dg  g(r,s,C)|  <  M  for  all  r  >  r0,  s  6  K,  E,  G  K'  and  |6|  <  m  +  5  (A. 10) 

Then  there  exists  a  constant  C  =  C(K,K',r0>g)  >  0  such  that  if  q(r,5)  is 
defined  by 


I(r,C)  -  (2-n)m'2 


exp{ i  r  0(t(5)  ,€)  +  sgn  9"(x(C)  ,5)  }g(r,x(0  .O 

r^2  | det  9"(x (0,5)  | 1/2 


+  q(r,C) 
(A. 11) 


then 


|q(r,C) |  <  C  r  ^  2  1  for  all  r  >  r0  and  £  e  K' 


(A. 12) 


In  (A. 11),  sgn  9"(s,£)  is  the  signature  of  the  real  symmetric 
matrix  0"(s,?). 

The  uniform  estimates  given  above  are  due,  in  essence,  to 
M.  Matsumura  [6].  No  proofs  are  offered  because  the  theorems  can  be 
proved  by  following  Matsumura' s  proofs  and  recognizing  that  under  the 
hypotheses  formulated  above  his  estimates  are  uniform  for  £  €  K'. 
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